Multi-magnon bound states in the frustrated ferromagnetic ID chain 
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We study a one-dimensional Heisenberg chain with competing ferromagnetic nearest-neighbor 
and antiferromagnetic next-nearest neighbor interactions in magnetic field. Starting from the fully 
polarized high-field state, we calculate the dispersions of the lowest-lying n-magnon excitations and 
the saturation field (n = 2, 3, 4). We show that the lowest-lying excitations are always bound multi- 
magnon states with a total momentum of 7r except for a small parameter range. We argue that bose 
condensation of the bound n magnons leads to novel Tomonaga-Luttinger liquids with multi-polar 
correlations; nematic and triatic ordered liquids correspond to n — 2 and n = 3. 

PACS numbers: 75.10.Jm, 75.10.Pq 



Frustrated spin chains are simple models that still show 
surprisingly rich physics. A prototype of these is the spin- 
i Heisenberg chain with nearest-neighbor (NN) J\ and 
next-nearest-neighbor (NNN) J2 couplings, 



H — ^ (^JiSi ■ Si + i + J2S1 ■ Si + 2 — hS'i 



(1) 



where Si = (Sf,Sf,Sf) is the spin- 1/2 operator on the 
ith site, and external field h is applied in the z direc- 
tion. The exchange interactions are frustrated as the 
NNN interaction is antiferromagnetic (AF), J2 > 0. Un- 
til recently most theoretical studies on the J\-Ji model 
(JTJ) considered the case where the NN coupling J% is also 
AF. However, interest is now growing rapidly in the fer- 
romagnetic (FM) case (Ji < 0) as well, which is triggered 
by experimental reports on thermodynamic properties of 
various quasi-one-dimensional frustrated FM spin chains 
(for a list of frustrated quasi-onc-dimcnsional materials, 
see Table 1 in Ref. [l| and Fig. 5 in Ref. |2|). For ex- 
ample, Rb 2 Cu 2 Mo 3 Oi2 (Ref. Q) and LiCuV0 4 (Ref. S) 
are considered to be described by the J1-J2 model with 
Ji w — 3J2 and Ji w — O.3J2, respectively. 

Recent theoretical studies 4 ^ have shown that, in a 
magnetic field, the FM Ji-J 2 chain (J 1 < 0, J 2 > 0) 
becomes a Tomonaga-Luttinger (TL) liquid having ne- 
matic quasi-long-range order as dominant correlation for 
some range of J\j Ji- This nematic state can be thought 
of as arising from condensation of two-magnon bound 
states^ Interestingly, such a nematic ordered phase can 
also appear in a frustrated FM spin model on the two- 
dimensional square latticed One can imagine further a 
bose condensed phase of more-than-two-magnon bound 
states. Indeed it was shown 8 - recently that an octupolar- 
like triatic ordered phase can result from condensation 
of three-magnon bound states in a frustrated ferromag- 
net on the triangular lattice. These results motivated us 
to examine the possibility of such many-magnon bound 
states in the FM J1-J2 model (JTJ. In this paper we show, 
by explicitly constructing bound-state wave functions, 
that for —3.52 < J1/J2 < —2.72 the lowest-lying exci- 
tations from the fully polarized FM state are 3-magnon 
bound states, and moreover 4-magnon bound states ap- 
pear for a slightly stronger FM coupling regime. We 



then suggest exotic TL liquid phases with multipolar-like 
spin correlations to emerge from condensation of multi- 
magnon bound states below the saturation field. 

Before presenting our calculations, let us briefly review 
known results that are relevant to our study. First, in the 
classical limit, we may regard Si as a c- number vector 
of length S. The ground state of (JJJ in this limit has 
a (right- or left-winding) helical spin structure, Si/S = 
(sin 8 cos (j>i, sin 8 sin <j>i, cos 6), with a pitch angle of <f> = 
<f>l+i — 4>i = ± arccos (— Ji/4J 2 ) and a canting angle of 
9 = arccos [Ah J2/ S(Ji + 4J2) 2 ], so that spins are fully 
polarized at a saturation field of h s = S(Ji + AJ2) 2 /4J2 
if -4J 2 < Ji < 0, or at zero field if J\ < -4J 2 . The 
helical order does not survive quantum fluctuations. 

In the quantum spin-i case, the ground state is fully 
polarized without magnetic field if J\ < — 4J2. At the 
boundary J\ = — 4 J2 the zero-field ground state is highly 
degenerate such that states from vanishing magnetiza- 
tion to full polarization share the same energy^ For 
the parameter range — 4J2 < J\ < of our main in- 
terest, Chubukov^ suggested that the ground state just 
below the saturation field should be a nematic state 
made up of bound magnon pairs with a commensu- 
rate total momentum k — tt if — 2.67J 2 < J\ < 
(2.67 w 1/0.38) and with an incommensurate momentum 
k < 7r otherwise, which was partly verified by mean-field 
theory,- numerical study^ Green's function analysis 
which fixed the commensurate-incommensurate transi- 
tion point to J1/J2 = -2.66908 (= -1/0.374661)^ and 
weak-coupling bosonization analysis^ ' 11 ' 13 While earlier 
calculations of the ground-state magnetization process 
suggested metamagnetic transitions ) 10 ! 14 recent density- 
matrix renormalization group (DMRG) study^i finds 
that the total magnetization of finite-size chains changes 
by AS Z = 2 at Ji = -J 2 , AS Z = 3 at Ji = -3J 2 , and 
AS Z = 4 at Ji = — 3.75J2 below saturation, implying 
that the magnetization curve is continuous in the ther- 
modynamic limit i^ 3 - 

To reveal how the fully polarized FM state collapses 
into new states with decreasing either magnetic field or 
the coupling ratio \J\\/ J2, we analyze magnon instabil- 
ity in the fully polarized state. We apply a large enough 
magnetic field h such that the fully polarized state is 
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the unique ground state for — 4J 2 < J\ < 0, and all 
multi-magnon excitations have positive excitation ener- 
gies which decrease as h is reduced. The saturation field 
h s is then defined as the field at which the lowest excita- 
tion energy vanishes. We consider bound states of up to 
four magnons in a finite chain of N spins. 

To explain our computational scheme, we begin with 
one- and two-magnon states. The one-magnon excited 
state with momentum k, 

1 N 

l fc ) = 7^E e ^ 5 r|FM), (2) 
V 1=1 

on the fully polarized state |FM) = | T T T - - - T ) l ias the 
excitation energy 

ei (k) = Ji(cosfc - 1) + J 2 [cos(2/c) - 1] + h, (3) 

which has a minimum of 



ei(fco) = - 




at cos fco = — </i/4J 2 . Obviously the one-magnon insta- 
bility just reproduces the classical helical spin state in 
the applied field below the saturation field. However, as 
shown by earlier studies;^ ' 10 ! 11 ' 12 this is not the true 
instability for the quantum case in the whole parameter 
region — 4J 2 < J\ < 0. 

To calculate bound n-magnon excitations, we take n- 
magnon states with a center-of-mass (CM) momentum k 
as a basis. For example, our basis for the two-magnon 
excitations has two [ spins with a total momentum k and 
a relative distance r (=1,2,.. .), 

|r; k) = -L J2 e lk{2l+r)/2 S-S- +r |FM) . (5) 

In this basis the matrix elements of the Hamiltonian 
can be written as (r; k\H\r' ; k') — Sk,k>H rtr r, where non- 
vanishing entries of H r ^ r i are 

if r r = Ji{8 r s — 2) + J 2 (<5r,l cos k + 8 r ^ — 2) + 2h, 
H r , r +i = H r+ i. r = Ji cos(fc/2), (6) 

H r ,r+2 — H r +2,r — ^2 COS/c. 

By separating off the CM motion we have reduced the 
two-magnon problem to a one-particle one which in prin- 
ciple can be solved exactly. For general k this involves 
finding roots of a transcendental equation.— The eigen- 
value problem ([6]) is greatly simplified at k = tt, where the 
excitation energy of the two-magnon bound state is&i 2 - 

72 

e 2 (7r) = -Ji-3J 2 + _ 2 _ + 2h. (7) 

Jl — J 2 

To calculate the energy dispersion e 2 (&) of the two- 
magnon bound states for general k S [0, 7r], we numeri- 
cally diagonalize the matrix H r y by restricting r and r' 



up to 1000. As long as the maximum value of r is suf- 
ficiently larger than the size of bound states, finite-size 
corrections should be exponentially small. 

The bound states of more than two magnons can be 
calculated in a similar manner. For the n-magnon sector, 
we take, as a basis set, states with total momentum k in 
which n down spins are separated by distance r±, . . ., 
r n _i. For example, the 3-magnon basis is given by 

ik{3l+2 ri +r 2 )/3 

l r *' r * fc > = E ^ S-S- +ri S- +ri+r2 |FM) , 

(8) 

with which the matrix elements (r±, r 2 ; k\H |r^, r 2 ; k) are 
easily found. The 4-magnon basis states \ri,r2,Ts]k) 
are constructed similarly. To solve the 3- and 4-magnon 
bound states, we numerically diagonalized the Hamilto- 
nian matrix expressed in terms of the finite number of 
basis states |ri,r 2 ;fc) with 1 < r-j < 27 and \ri, r 2 , r 3 ; k) 
with 1 < Ti < 9, respectively. This seems to be sufficient 
to determine the lowest excitations. 

The four panels of Fig. Q] show energy dispersions of 
multi-magnon excitations at Ji/J 2 = —1.0, —2.7, —3.0, 
and —3.8. For each value of J\/ J% the magnetic field is 
set equal to the saturation field where the lowest-lying ex- 
cited state is gapless. The dispersion of bound n — 2, 3, 4 
magnons are obtained by diagonalizing the Hamiltonian 
in the finite basis, as we described above. Our numeri- 
cal calculation reproduces the known results for the two- 
magnon sectori 6 i 10 i 12 the lowest bound two-magnon exci- 
tation has momentum k = tt for — 2.67 < Ji/J 2 < and 
incommensurate momentum k < tt for —4 < Ji/J 2 < 
—2.67. Figure Q] also shows the lower edges (thin lines) 
of the continuum spectra of scattering states made up 
of n magnons, which are calculated from the dispersion 
of a single magnon §5§ and those of bound states of up 
to n — 1 magnons. For example, energies of 3-magnon 
scattering states are obtained by adding those of either 
three unbound magnons or of a 2-magnon bound state 
plus a free magnon, under the condition that the total 
momentum is k (mod 2tt). 

We see in Fig. [TJthat for any J\l J% there is always a 
region in the k space where bound states lie well below 
the scattering continuum. The character of the lowest- 
lying bound states signalling the instability of the fully 
polarized state changes with Ji/J 2 . Unlike previously 
thought, the system shows rather different regimes as a 
function of Ji/J 2 : a two-magnon commensurate (k = tt) 
instability (—2.67 < Ji/J 2 < 0), a two-magnon incom- 
mensurate instability as predicted in Ref. [6; (—2.72 < 
J\/ J2 < —2.67), a three-magnon commensurate instabil- 
ity (—3.52 < J1/J2 < —2.72), and a four-magnon com- 
mensurate instability (—3.86 < Ji/J 2 < — 3.52)»^ As 
Ji/J 2 approaches further towards —4, we expect to have 
more-than- four-magnon commensurate phases, but these 
are outside the scope of our numerics. We note that, ex- 
cept in the incommensurate regime (—2.72 < J1/J2 < 
—2.67), the bottom of the lowest-lying mode is at k = tt. 
For example, at J1/J2 — —1.0 the lowest energy state 



3 





; Ji/J 2 =-3.0 



.y<- 







0.25 



0.5 

k/K 



0.75 




0.075 - 



0.05 



0.025 



FIG. 1: (Color online) Energy dispersions for the multi-magnon bands at the saturation field. Thin lines denote the onset 
of the scattering continuum, and thick lines show n-magnon bound states where relevant. With decreasing J1/J2, the lowest 
excitations change from commensurate (k — n) two-magnon bound states, incommensurate (k < n) two-magnon bound states, 
commensurate three-magnon bound states, to commensurate four-magnon bound states. 



in the 2-magnon sector is the bound state with k = tt. 
The 4-magnon scattering states of two such two-magnon 
bound states therefore have lowest energy (which vanish 
when h = h s ) at k = 0. The fact that we do not have 
a mode of 4-magnon bound states near k = below the 
continuum indicates that the interaction between two 2- 
magnon bound states of k — tt is repulsive, and therefore 
these 2-magnon bound states are stable. At J\ = — 3J 2 
we have the 3-magnon bound state with k = tt as the 
lowest-energy excitation. The presence of stable 3- and 
4-magnon bound states provides natural explanation for 
the AS Z > 2 jumps found in the DMRG study of Ref.Qj]. 

Figure [5] shows the saturation field h s determined by 
the instability from the softening of the lowest exci- 
tations. The saturation field estimated from a single- 
magnon instability is always smaller than the true satu- 
ration field which is determined by multi-magnon bound 
states. The calculated saturation field is in perfect 
agreement with the exact h s estimated from Eq. ([7]) for 
-2.67 < J1/J2 < 0. For -3.86 < J1/J2 < -2.72 the 
saturation field is determined from the instability of 3- 
or more-magnon bound states. 



In Fig. [3] we show the expectation value (^27=1 Ti ) 
characterizing the size of the n-magnon bound states, 
where the average is taken for the lowest-energy bound 
state of n magnons for a given Ji/ J 2 ; the minimum value 
of the average is by definition n — 1 for the n-magnon 
bound state. We confirm that the magnons are tightly 
bound when they are the lowest-lying excitations in the 
energy spectra. This justifies our use of basis states with 
finite Ti for calculating the low-lying bound states. 

We now discuss the implications of the multi-magnon 
instabilities that lead to condensation of bound magnons 
below the saturation field. We argue that the bound 
n-magnon condensation gives rise to a phase with 
multipolar-like quasi-long-range order but without spin 
(dipole) ordering in the XY direction. To be specific, let 
us consider the case where the n-magnon bound states 
with momentum k = tt become gapless as h — > h s + 0; 
n = 2 for -2.67 < J1/J2 < 0, n = 3 for -3.52 < Ji/J 2 < 
-2.72, and n = 4 for -3.86 < Ji/J 2 < -3.52. Just be- 
low the saturation field, the system can be viewed as a di- 
lute gas of repulsively interacting bosons which represent 
the n-magnon bound states^ Here the boson creation 
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FIG. 2: (Color online) Saturation field versus J1/J2. The 
solid curve "n magnon" shows the saturation field h s obtained 
from the numerical n- magnon solutions (n = 2,3,4), while 
the curves "1 magnon" and "2 magnon" denote h s calculated 
from Eqs. ((4| and respectively. The inset shows the re- 
gion where our saturation field deviates from the two-magnon 
solution and the saturation fields for each multi-magnon ex- 
citation. 
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FIG. 3: (Color online) Mean length r^} of the bound 
multi-magnon states versus Ji / J2 . The kinks in the curves ap- 
pear when the total momentum of the n-magnon bound state 
changes between incommensurate (where (Y]- rj) is larger) 
and commensurate (k — tt). 



operator bj is identified with (— lySfS^ 
crude approximation. These bosons condense to form a 
TL liquid^ 6 - in which various correlation functions of the 
boson fields decay algebraically. The most slowly decay- 
ing two-point correlation function will be the propagator 
of the bosons (bob}.), which decays as (— l) r r _1 / ?? with 17 
77 — ► 2 as h — ► h g — 0. However, the spin operators 
which cannot be simply represented with bosons bi, such 
as Sf and its products IC^SJ" with p = 0, 1, • • ■ , n — 2, 
should have only short-range correlations. This is be- 
cause exciting one magnon costs binding energy. Since n 
down spins form a tightly bound state, we may use the 
approximation A — Sf — nb\b i . This allows us to cal- 
culate the correlation function of Sf from the density 
correlation of bosons, which has the asymptotic form 
{b\b blb r ) ~ p 2 + Ar~ n cos(27rpr) — r/(2irr)~ 2 , where p 
is the boson density (np = \ — (S z )) and A is a con- 
stant. For the case n = 2, the above theory indicates 
that the TL liquid has nematic quasi-long-range order, 6 
which is indeed found by the recent DMRG calculation at 
J 1 = — J2.— ^ The theory also predicts that the TL liquid 
with larger n (n = 3,4) should exist for larger \J\\/ J2 in 
the phase diagram of the FM J\- J2 model near the satu- 
ration field. For n — 3 this is a TL liquid with antiferro- 
triatic quasi-long-range order, which is a one-dimensional 
analogue of triatic order found in Ref. Our numerics 
suggests that, as J1/J2 approaches —4, instability from 
bound states of more magnons appears. We do not know 
how far these new phases extend to lower fields. 



In summary we have numerically calculated many- 
magnon bound states and determined their energy dis- 
persions. We have found that the fully polarized FM 
state has instabilities to bose condensation of these many- 
magnon bound states, which lead to TL liquids with mul- 
tipolar magnetic correlations below the saturation field. 
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The value of J\ = — 3.86J2 for the lower edge of the 4- 
magnon region is a lower bound. At Ji < — 3.86J2 a loosely 
bound incommensurate 4-magnon state has lower energy 



than the commensurate bound state. It is quite possible 
that 5-magnon bound states with k — tt becomes lowest 
excitation at a slightly larger value, J1/J2 ~ —3.8. 
Assumption of mutual repulsive interaction is justified by 
the absence of a metamagnetic transition below the satu- 
ration field and by the AS = n jumps in the magnetization 
curve . 11 ' 13 For the repulsive interaction the TL parameter 
K is less than 1. 

The exponent 77 is related to K as r\ = 2K. Reference 
obtained 77 = 1.1 for Ji = -J 2 (n = 2) and (S z ) = 3/8. 



